Chandra detection of diffuse hot gas in and around the M31 bulge 
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. We report the detection of diffuse hot gas in M31, using archival Chandra ob- 

servations which allow us to map out a 30' x 30' field (covering a galactocentric 

r^ i radius up to 4.5 kpc) and to detect sources in the galaxy down to a 0.5-8 keV lu- 

minosity limit of ~ 10 35 ergs s _1 . We estimate the remaining stellar contribution 
from fainter X-ray sources (primarily cataclysmic variables and coronally active 

£j [ binaries), assuming that they spatially follow the stellar distribution. Indeed, the 

near-IR K-band light of the galaxy closely traces the 2-8 keV unresolved X-rays, 
indicating a collective stellar X-ray emissivity consistent with those determined 
for the Galactic ridge and M32, whereas the amount of the 0.5-2 keV unresolved 
emission is significantly greater than the expected stellar contribution, especially 

[ — , 

within a galactocentric radius of ~2 kpc. Morphologically, this soft X-ray excess 
appears substantially rounder than the bulge as seen in K-band and is elongated 
approximately along the minor-axis at large radii. The excess thus most likely 

t — j represents the emission of diffuse hot gas in and around the galactic bulge. Fur- 

thermore, the near side of the M31 disk casts an apparent shadow against the 

.{ph ! soft X-ray excess, indicating that the hot gas extends to at least 2.5 kpc from 

the galactic plane. We briefly discuss the implications of these results on the 
energy balance in the M31 bulge and on understanding the large-scale soft X-ray 
enhancement observed toward the inner region of our own Galaxy. 

Subject headings: galaxies: general — galaxies: individual (M31) - galaxies: 
spiral — X-rays: general 



1. Introduction 

The bulge of a galaxy, though consisting of mainly old stars, is a mecca of high-energy ac- 
tivities. Low-mass X-ray binaries (LMXBs), in the luminosity range of ~ 10 35 — 10 38 ergs s^ 1 , 
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are among the brightest X-ray sources observed in a galaxy like our own. At lower lumi- 
nosities, typically ~ 10 30 — 10 33 ergs s" 1 , are numorous cataclysmic variables (CVs) and 
coronally active binaries (ABs) that can be individually detected normally only in the Solar 
neighborhood. Such stars are most likely responsible for the bulk of the unresolved 2-10 keV 
emission observed in the Galactic bulge/ridge (Revnivtsev et al. 2006), whereas lower energy 
X-rays from the same regions are subject to heavy interstellar absorption and hence diffi- 
cult to detect. One thus needs external perspectives of nearby galaxies. Indeed, Revnivtsev 
et al. (2007) have shown that the unresolved emission from the low-mass bulge-dominated 
galaxy M32 over the entire 0.5-7 keV range is primarily stellar in origin. 

The interstellar medium (ISM) in a galactic bulge is also expected to be extremely 
energetic, chiefly due to a concentration of Type la supernovae (SNe). The bulk of the 
mechanical energy release from such SNe is expected to be in shock-heated gas which can be 
naturally traced by its X-ray emission. However, the observed luminosity of the unresolved 
(source-removed) X-ray emission from a galactic bulge typically accounts for only a small 
fraction (a few %) of the expected SNe energy release. How the remaining energy is dissipated 
remains unknown. It may be propagated into the large-scale halo of the host galaxy in a 
mechanical outflow or sound waves, for example. In any case, determining the fate of this 
"missing" stellar feedback energy is fundamentally important in our understanding of its 
role in galaxy evolution (e.g., Wang 2007). Here we present the first step of such a study. 

We detect the truly diffuse soft X-ray emission from the bulge of M3 1 — the nearest spi- 
ral galaxy (d ~780 kpc; l'=0.23 kpc) that is similar to our own Galaxy. The galaxy contains 
no AGN, so the non-nuclear X-ray emission can be studied even in the very central region. 
The M31 bulge has little cool gas and star formation, so the X-ray contribution from young 
stellar populations is minimal. The moderate inclination (~ 78°) of the M31disk and the 
relatively low Galactic foreground absorption (N H ~ 6.7xl0 20 cm' 2 ) also allow us to detect 
extraplanar X-ray emission in the 0.5-2 keV range. Indeed, detections of diffuse soft X-ray 
emission have been claimed, based on spectral decompositions (Shirey et al. 2001; Takahashi 
et al. 2004). However, such a decomposition depends sensitively on rather arbitrary choices 
of spectral models for various components contributing to the spectrum, both diffuse and 
discrete. Our approach here is to spatially map out the diffuse X-ray emission (in addition to 
its energy dependence) and to study its relationship to other stellar and interstellar compo- 
nents of the galaxy. While details of the data reduction as well as a comprehensive analysis 
of the sources and modeling of the hot gas will be presented elsewhere (Li et al. 2007, in 
preparation), we focus here in presenting the detection of the diffuse hot gas in and around 
the M31 bulge and its immediate implications. We quote our statistical error estimates at 
the 90% confidence level. 
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2. Data preparation 

Our X-ray study was based on 31 Chandra ACIS archival observations of M31 taken 
by 2005. The majority (21 out of 31) of these observations were taken with the ACIS- 
I array and aimed toward the M31 bulge with the aim-points located within 1' from the 
galactic center. To maximize the coverage and uniformity of the combined field, we utilized 
data only from the front-illuminated CCDs (the ACIS-I array and the S2 chips) of the 21 
observations. For same reason, we also included I-chip data from four ACIS-S observations. 
These data together cover a field of r ~18' around the center of M31. Furthermore, for 
local sky background determination, we used six additional ACIS-I observations which were 
aimed toward an "off-field" ~20' southwest to the center. 

We reprocessed the data using CIAO (version 3.3), following the Chandra ACIS data 
analysis guide. We generated count and exposure maps for each observation in the 0.5-1, 1-2, 
2-4 and 4-8 keV bands. Corresponding instrumental background maps were generated from 
the "stowed" data, after calibrating the 10-12 keV count rate with individual observations. 
The total effective exposure is ~95 ks in the central region and gradually drops to < 20 ks 
at radii r > 10'. 

Following a procedure detailed in Wang (2004), we performed source detection in the 
soft (0.5-2 keV), hard (2-8 keV) and broad (0.5-8 keV) bands. With a local false detection 
probability P < 10~ 6 , a total of 305 sources are detected in the field. To study the unresolved 
X-ray emission, we excluded each of the sources from maps of individual observations with 
circular regions enclosing ~97% of the source counts. The residual of this source removal 
contributes about 10% of the remaining unresolved X-ray emission in the field. The source- 
removed maps were then reprojected to generate combined images in the four bands. We 
further statistically corrected for the variation of the detection incompleteness across the 
field, to a common detection limit of 8xl0 34 ergs s _1 (0.5-8 keV). Because of the relatively 
flat luminosity function of the sources (mostly LMXBs; Li et al. 2007 in preparation; see 
also Voss & Gilfanov 2007), the correction (normalized according to the 2MASS K-band 
intensity; Fig. [Th.; Jarrett et al. 2003) typically amounts to less than 6% of the unresolved 
emission. For the same reason, the residual contribution from LMXBs at lower luminosities 
is found to be negligible. 

3. Analysis and results 

Fig. [T^i shows the 0.5-2 keV unresolved X-ray emission from a 30' by 30' region around 
the center of M31, compared with the near-IR K-band image. In the inner bulge and along 
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the major-axis of the disk, the X-ray emission shows morphological similarities with the ri- 
band light, whereas at large radii the X-ray morphology appears substantially rounder and is 
elongated approximately along the minor-axis, in particular at the southeast side, indicating 
the presence of diffuse hot gas. 

3.1. The collective stellar emission 

To quantify the diffuse hot gas, we need to isolate the collective stellar contribution, 
which presumably spatially follows the K-band light distribution. Optimal for this purpose 
is to inspect regions along the major-axis, where the soft X-ray emission morphologically 
mimics the K-band light better than in regions further away from the major-axis (Fig. [T^l) . 
Fig. [2^ shows the unresolved X-ray intensity profiles along the major-axis, together with 
the corresponding K-band intensity profile. Indeed, the K-band profile as a "model" fits 
the hard band profile (triangles in Fig. [2h) well with a normalization factor Nk = 4.3 ± 
0.2 xl0~ 5 cts s _1 arcmin~ 2 /(MJy sr _1 ). Therefore, the hard band X-ray emission is fully 
consistent with an origin in the old stellar population. 

The collective stellar emission should also contribute at lower energies. However, while 
the soft X-ray profile (diamonds in Fig. [2h) can match that of the K-band light reasonably 
well at major-axis radii > 8', there is a clear excess above the collective stellar contribution 
in the inner region. This soft excess in the bulge is another indication for the presence 
of hot gas, although its exact spatial distribution is yet to be determined. We include an 
exponential law to approximately account for the excess. The resultant fit is satisfactory 
(solid curve in Fig. Hi), with a fitted N K = (25.0 ± 2.4) x cts s" 1 arcmin" 2 /(MJy sr" 1 ). 
The fit predicts that the stellar component contributes about 60% to the soft emission within 
a major-axis raidus of ~8' and becomes dominant further beyond. The fit also reveals an 
interesting drop of the X-ray to K-band intensity ratio within the central 0'5 (more clearly 
indicated in a radial intensity profile which is not shown here) , with values still comparable to 
or greater than that at the large major-axis radii. The nature of this ratio drop is currrently 
being investigated, and we note that it has little effect on the result of the above procedure. 

3.2. The diffuse X-ray emission 

Fig. [lb shows an image of the truly diffuse emission after subtraction of the collective 
stellar contribution from the total unresolved X-ray emission of M31. The emission along the 
major-axis is confined within a projected distance of ~8' (~1.8 kpc) from the galactic center 



- 5 - 



to the southwest and is slightly more extended to the northeast. The overall morphology is 
elongated approximately along the minor-axis, with an extent of more than 15' (~3.5 kpc) 
on both sides with respect to the center; but the emission appears considerably fainter on the 
northwest side and somewhat interrupted by the presence of spiral arms. This asymmetry 
is further illustrated in Fig. [2b, where diffuse X-ray intensity profiles along the minor-axis 
are shown separately in the 0.5-1 keV and 1-2 keV bands. The asymmetry can be naturally 
explained by the soft X-ray absorption of the galactic disk, as its near side is to the northwest. 
In particular, a major spiral arm and the star-forming ring, as traced by the peaks of the 
Spitzer MIPS 24 /zm emission (Gordan et al. 2006), apparently cast deep X-ray shadows on 
the northwestern side. Estimated from the relative depth of these shadows, the equivalent 
X-ray- absorbing column densities are ~1.2xl0 21 cm -2 and ~3.6xl0 21 cm -2 , consistent with 
the hydrogen column densities of the spiral arm and the star- forming ring (Nieten et al. 2006). 
No similar shadow is apparent on the southeastern side, indicating that the emission on this 
side is mostly from the bulge region in front of the disk. Therefore, the diffuse emission 
seems to have an intrinsic (absorption-corrected) overall coherent morphology reminiscent 
of a bi-polar outflow from the bulge. The soft X-ray absorption by the disk suggests that 
the vertical extent of the X-ray-emitting gas from the galactic plane is at least 2.5 kpc. 

We characterize the profiles (Fig. [2b) at distances of —7' < z < 0' off the major-axis 
with an exponential law: I(z) = I g e~^ z ^ z °, where I g is the central intensity and zq is the 
projected scale-height. The best-fit z , being 2f5=b0f 1 (~0.6 kpc), shows no statistically 
significant difference between the two bands, indicating that the hot gas in the bulge has 
little temperature variation. The hardness ratio, / gi i_2fc e v/-^,o.5-ifceV~0.25, is consistent 
with a spectrum from an isothermal gas with a temperature of ~0.4 keV, subject to the 
Galactic foreground absorption. For comparison, the hardness ratio of the stellar emission 
is ~0.72; the scale-heigh of the K-band light is ~2'1, if characterized by an exponential 
law as well. Therefore, the diffuse X-ray emission is both softer and more extended than 
the stellar contribution. At \z\ > 7', however, the soft X-ray intensity distribution levels 
off (Fig. [2b). Part of this leveling may be related to emission associated with the galactic 
disk, partially compensating its X-ray absorption. The exact intensity level at such large 
distances, however, depends on an accurate subtraction of the local sky background which 
has been estimated in the off- field and may be biased (§ 2). The sky coverage of the present 
ACIS observations is still too limited to accurately determine both the background and the 
large-scale distribution of the diffuse emission associated with the bulge. Assuming the above 
exponential fit and temperature estimate as well as an intrinsic symmetry with respect to the 
galactic plane, we infer a 0.5-2 keV luminosity of the diffuse emission as ~2.2x 10 38 ergs s -1 . 
The flat tail parts of the profiles give an additional ~3xl0 37 ergs s" 1 . 
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4. Discussion 

In § !3.1l we have estimated the stellar contribution to the total unresolved X-ray emission, 
which is a key step in isolating the truly diffuse emission. It is thus instructive to compare 
our result with independent measurements. Sazonov et al. (2006) measured the collective X- 
ray emissivity (per unit stellar mass) of the old stellar populations in the Solar neighborhood 
to be 9±3 x 10 27 ergs s" 1 M" 1 in the 0.5-2 keV band and 3.1±0.8 x 10 27 ergs s" 1 M" 1 in 
the 2-10 keV band. Revnivtsev et al. (2006) showed that the Galactic ridge X-ray emission 
closely follows the near-IR light that traces the Galactic stellar mass distribution, and that 
the X-ray to near-IR intensity ratio is consistent with the collective X-ray emissivity of old 
stellar populations inferred from the Solar neighborhood. Revnivtsev et al. (2007) further 
found that the 0.5-7 keV unresolved X-ray emission and K-band stellar light in M32 have 
consistent spatial distributions, but they did not explicitly give fitted parameters of the 
spectral model, which would allow for an immediate comparison with the M31 values. 

We have thus re-extracted the unresolved X-ray spectrum of M32 from two Chandra 
ACIS-S observations (Obs.ID. 2017 and 5690) with a total exposure of 160 ks. The spec- 
trum can be adequately fitted by a model consisting of a power-law component (with a 
photon index of 1.86^q|i) and a thermal plasma component (temperature of 0.45^o;Jo keV), 
with the Galactic foreground absorption. The collective emissivity is 5.8±1.1(5.6±1.0) x 
10 27 ergs s -1 Mq 1 in the 0.5-2 (2-10) keV band, consistent with the values reported by 
Revnivtsev et al. (2007). The spectral model-predicted ACIS-I count rate is ~17. 1 (3.9) 
xlO -5 cts s~ 1 arcmin _2 /(MJy sr _1 ) in the 0.5-2 (2-8) keV band. The hard band value agrees 
well with the measurement for M31 (4.3±0.2x 10~ 5 ), accounting for the residual photons 
spilling outside the source- removal regions (~ 10%; § [2]). The soft band stellar emissivity 
of M31 (residual LMXB contribution excluded) is a factor of 1.4 ± 0.3 higher than the M32 
value, but is consistent with that inferred from the Solar neighborhood. This discrepancy 
is not totally unexpected, considering various statistical errors and hiden systematic uncer- 
tainties (e.g., in the mass-to-light ratio and in the spatial and spectral modeling). We note 
that adopting the M32 stellar emissivity to remove the stellar contribution in M31 would 
enhance the hot gas contribution in and around the M31 bulge, but would not qualitatively 
alter the picture of the diffuse emission presented in § 13.21 and below. 

The characterization of the diffuse hot gas sheds important insights into the energy 
balance in the M31 bulge. The estimated luminosity of the hot gas (2.5 x 10 38 ergs s _1 ) 
is only about 0.6% of the expected SNe mechanical energy input, ~4 x 10 40 ergs s _1 . As 
mentioned in § [TJ this indicates that the input energy may be removed primarily in an 
outflow. Dynamically, such an outflow tends to find its way along steeper pressure gradient 
against the gravity of the galaxy, consistent with the observed bi-polar morphology of the 
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diffuse X-ray emission. If the gas were quasi-static, one would expect its distribution to follow 
that of the gravitational potential, i.e., more extended along the major-axis. However, the 
gas may not be hot enough to ultimately escape from the deep gravitational potential of 
M31; it is also not clear how the outflow interacts with the large-scale halo of M31 and 
how the mechanical energy is dissipated. Similar considerations also challenge the studies of 
X-ray- faint elliptical galaxies (e.g., David et al. 2006). Ongoing numerical simulations would 
help to understand the nature of the hot gas and its role in the evolution of these systems. 

Our unambiguous detection of the diffuse hot gas in and around the M31 bulge also 
helps to understand the soft X-ray enhancement observed toward the inner region of our 
Galaxy. The temperature of the hot gas associated with the M31 bulge, 0.4 keV, is similar 
to that with the Galactic bulge, as estimated from the ROSAT all-sky survey (Snowden 
et al. 1997). Based on a hydrostatic model of the Galactic bulge X-ray emission devel- 
oped by Wang (1997), Almy et al. (2000) further inferred a total 0.5-2 keV luminosity of 
~ 8 x 10 38 ergs s _1 , about four times greater than our estimated M31 bulge luminosity. 
The relatively high luminosity of the Galactic bulge manifests in the large extent of the 
soft X-ray enhancement from the Galactic bulge. At Galactic lattitudes b ~ —15° (~2 
kpc from the plane), for example, where both the confusion with the foreground emission 
features and the interstellar absorption are relatively small, the intensity has an averaged 
value of ~6(4)xl0 -4 ROSAT PSPC cts s -1 arcmin" 2 in the 0.75 (1.5) keV band (Snowden 
et al. 1997). Had this emission been detected from M31 by the Chandra ACIS-I, it would be 
measured with an intensity of ~10(3)xl0~ 4 cts s _1 arcmin -2 in the 0.5-1 (1-2) keV band, 
about 2-4 times higher than the observed M31 values represented by the tails (Fig. [2b). 
The intensity drops slowly and even shows local enhancements at high latitudes (b > —30°; 
Snowden et al. 1997). Within |6| < 10°, the interstellar absorption is severe, little can be 
inferred reliably about the properties of the hot gas. It is in this corresponding region in 
the M31 bulge (\z\ < 6') that the diffuse soft X-ray intensity shows the steepest increase (by 
about one order of magnitute) toward the galactic center. Such a mid-plane concentration of 
diffuse soft X-ray emission may also be present instrinsically in our Galactic bulge. Clealy, 
a more careful comparison and modeling of the X-ray data sets are needed in order to un- 
derstand the similarity and difference in the hot gas characteristics and their relationship 
to other galactic properties (e.g., the effect of recent active star formation in the Galactic 
center) . 

We thank D. Calzetti, M. Fardal, and S. Tang for helpful comments and discussions. 
This work is supported by the SAO grant AR7-8006X. 
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RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Fig. 1. — (a) Background-subtracted, exposure-corrected and smoothed intensity contours of 
the Chandra ACIS-I 0.5-2 keV unresolved emission overlaid on the 2MASS K-band image of 
M31. The contours are at 3, 6, 10, 16, 32, 64, 128, 196, 256 and 512 xl0~ 4 cts s" 1 arcmin" 2 . 
The galactic center is marked by a plus sign, (b) Contours of the diffuse (stellar contribution- 
subtracted) X-ray intensity (solid) and K-band light (dotted) overlaid on the Spitzer MIPS 
24 /zm image. 
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Vertical Distance (arcmin) Vertical Distance (arcmin) 

Fig. 2. — (a) 0.5-2 (diamonds) and 2-8 keV (triangles) keV intensity profiles of the unresolved 
X-ray emission along the major-axis. A position angle of 45° is adopted. The positive side 
is toward the southwest. The full width for averaging the intensity is 8'. Spatial binning is 
adaptively ajusted to achieve a signal-to-noise ratio better than 4, with a minimum step of 
6". The X-ray profiles are characterized by a normalized K-band intensity profile (dashed 
curves), and an additional exponential law for the soft band (solid curve), (b) 0.5-1 (crosses) 
and 1-2 keV (diamonds) intensity profiles of the diffuse emission along the minor- axis; the 
stellar contribution has been subtracted. The positive side is toward the northwest. The full 
width for averaging the intensity is 16'. The adaptive steps achieve a signal-to-noise ratio 
better than 3. The solid curves represent a fit with an exponential law. The corresponding 
24 /im intensity profile is shown by the dashed curve. The arrows mark the positions of the 
shadows casted by a spiral arm and the star-forming ring. 



